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Phase diagram in an external magnetic field beyond a mean-field approximation 
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The phase structure of the Polyakov loop-extended chiral quark-meson model is explored in 
a nonperturbative approach, beyond a mean-field approximation, in the presence of a magnetic 
field. We show that by including meson fluctuations one cannot resolve the qualitative discrepancy 
on the dependence of the crossover transition temperature in a non-zero magnetic field between 
effective model predictions and recent lattice results (arXiv:f if 1.4956). We compute the curvature 
of the crossover line in the T — \ib plane at a non-zero magnetic field and show that the curvature 
increases with increasing magnetic field. On the basis of QCD inequalities, we also argue that, at 
least in the large N c limit, a chiral critical end point and, consequently, a change from crossover 
to a first-order chiral phase transition are excluded at zero baryon chemical potential and non-zero 
magnetic field. 
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I. INTRODUCTION 

The phase diagram of strongly interacting matter in 
the presence of an external magnetic field has been been 
explored in model studies and in first principle Lattice 
Quantum Chromodynamics (LQCD) calculations. These 
investigations has an importance for physics of heavy-ion 
collisions at top RHIC and LHC energies, where accord- 
ing to the estimates (see Refs. p] and [2]), the magnitude 
of the magnetic field may reach extremely large values of 
eB ~ XQml. 

Model studies (see e.g. Refs. [SHZ|) have revealed a gen- 
eral structure of the phase diagram in the temperature- 
magnetic field (T — eB) plane. They showed that the 
chiral crossover temperature T pc increases with increas- 
ing magnetic field. This result has been, however, ob- 
tained within a mean-field approximation, i.e. meson 
fluctuations have not been taken into account. How- 
ever, charged pion degrees of freedom interacting with 
the magnetic field may drastically change properties of 
the crossover transition even at moderate values of B. 

Early LQCD calculations [5] confirmed the increase of 
T pc with eB. Unfortunately, results of Ref. [S] may not 
be final: the calculations were performed with standard 
staggered fermions, large lattice spacing and with some- 
what heavy pion masses. These deficiencies may diminish 
the role of charged pions. 

Recent LQCD studies [H] at physical pion mass, with 
improved staggered fermions and extrapolation to the 
continuum have shown completely opposite dependence 
of T pc on the magnetic field eB. It has been found that 
the transition temperature significantly decreases with 
increasing magnetic field. 

In this article, by analysing the chiral structure of 
the phase diagram of the Polyakov loop-extended quark- 
meson model beyond a mean-field approximation, we 
demonstrate, that meson fluctuations cannot explain 
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above mentioned qualitative discrepancy in the behavior 
of the phase transition temperature with the magnetic 
field. 

Some model calculations have also predicted that at 
large values of the magnetic field the strength of the tran- 
sition increases and eventually, at a critical end point, 
the chiral crossover turns into a first-order phase tran- 
sition |5j. Our model calculations with physical number 
of colors N c = 3 show no evidence in favour of possible 
chiral first-order phase transition in a wide range of the 
magnetic field. Based on the no-go theorem formulated 
for QCD in Ref. [TU], we will also argue that a chiral 
critical end point and a change from chiral crossover to 
a chiral first-order phase transition are forbidden at the 
leading order of the large N c expansion for any non-zero 
pion mass. 

In this paper, we do not consider the, so-called, no- 
sea approximation, the essence of which boils down to 
neglecting vacuum contributions in a thermodynamic po- 
tential. As was shown analytically in Ref. jTT] , the no-sea 
approximation results in a first-order chiral phase transi- 
tion at zero chemical potential in the chiral limit (already 
on the mean-field level). This is in direct contradiction 
with recent LQCD findings [12 . We refer to Refs. [3j 
and [H] for the reader interested in a comparison be- 
tween results of a properly renormalized theory and those 
obtained in the no-sea approximation at finite and zero 
magnetic field correspondingly. 

The structure of the paper is as follows. In the next 
section we discuss the functional renormalization group 
flow equations in the non-trivial magnetic background. 
In Section 3, we review the mean-field approximations for 
the Polyakov loop-extended Quark-Meson (PQM) model. 
Since the PQM model is renormalizable, we perform di- 
mensional regularization to subtract divergences arising 
from the vacuum term. Section 4 is devoted to the main 
results and discussions, including above mentioned input 
from the large N c limit. Section 5 contains our conclu- 
sions. 
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II. THE POLYAKOV LOOP-EXTENDED 
QUARK-MESON MODEL 

The quark-meson model is an effective realization of 
the low-energy sector of QCD. The model is built to re- 
spect symmetries of QCD, including the chiral symmetry 
in the limit of vanishing pion mass. However, because the 
local color SU(N C ) invariance of QCD is replaced by a 
global symmetry, the model does not describe confine- 
ment. The improved version of the model, the so-called 
Polyakov loop-extended quark-meson (PQM) model in- 
corporates a coupling of the quarks to a uniform tempo- 
ral color gauge field, represented by the Polyakov loop. 
In the PQM model, many facets of confinement can be 
emulated [T3HT5] . 

The Lagrangian of the PQM model reads 

C = q [i-fDp - g(a + i 75 T7f)] q - *•) (1) 

The coupling between the effective gluon field and 
quarks, and between the (electro)magnetic field B and 
quarks is implemented through the covariant derivative 



(2) 



where = gA%\ a /2 and A™ = (0, Bx, 0, 0). The 
flavor matrix Q is defined by the quark electric charges 
Q = diag(2/3e, — l/3e). The spatial components of the 
gluon field are neglected, i.e. A^ = S^qAq. The inter- 
action of charged pion ^ — (m ± WT2)/-\/2 with the 
electromagnetic field is included by 2? M = <9 M + ieA® . 
The effective potential for the gluon field U(&, $*) is ex- 
pressed in terms of the thermal expectation values of the 
color trace of the Polyakov loop and its conjugate 



1 



(Tr c L(x)}, 



1 

N~< 



(Tr c it(f)), (3) 



with 



■Pexp 



drA^(S, t) 



(4) 



where V stands for the path ordering, j3 = 1/T and 
A 4 = iAq. In the 0(4) representation, the meson field 



is introduced as 4> — i* 7 ^ — (^o? tti, ^2)) and the corre- 
sponding SU(2)l ® SU(2)r chiral representation is de- 
fined by er + if ■ ir~f 5 . 

The meson potential of the model, U (a, ff) , is defined 

as 



A 



{/(cr,??) = - (<7 2 + n 2 -V 2 )' 



(5) 



The effective potential of the gluon field is parametrized 
in such a way as to preserve the Z(3) invariance: 



Z/($,$*) = b 2 {T) 
The parameters, 
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b 2 (T) = a Q + ai \ — j +a 2 
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T 



(6) 



(7) 



with a = 6.75, ai = -1.95, a 2 = 2.625, a 3 = -7.44, 
63 = 0.75, b 4 = 7.5 and T = 270 MeV were chosen to 
reproduce the equation of state of the pure SU C (3) lattice 
gauge theory. When the coupling to the quark degrees 
of freedom are neglected, the potential (JoJ) yields a first- 
order deconfinement phase transition at T - 



A. The FRG method in the PQM model 

To take fluctuations into account in the PQM model, 
we use the functional renormalization group (FRG) 
method. This method is based on an infrared regular- 
ization of the fluctuations at a sliding momentum scale 
fc, resulting in a scale-dependent effective action Tk, the 
so-called effective average action [16J . In this article, the 
Polyakov loop is treated as a background field introduced 
self-consistently on the mean-field level. Quark and me- 
son fluctuations are accounted for by solving the FRG 
flow equations. 

The FRG flow equation for the PQM model at zero 
magnetic field was derived in Ref. [17j . The derivation 
for finite magnetic field is lengthy, but straightforward 
and can be done along similar lines. Here we provide 
only the final expression for the flow equation. It can be 
easily proved that in the limit B — > the flow equation 
reduces to that of Ref. [T7] . 



Following our previous work |17) . the flow equation for the scale-dependent grand canonical potential density, 
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fife = TTk/V, is formulated for the quark and meson subsystems in the leading order of derivative expansion 



12tt 2 



2tt 2 



l + 2n B ( J E ;T ;r) 
l + 2n B (£ 7r ;T)' 



(8) 



E E v fc2 - «i( n > ( 3//' s ) 6l ( fc2 - 0//. *)) 



2tt 2 

/=1,2 s= l_i n =0 



Here nsiE^^; T) is the bosonic distribution function 
with the pion and sigma energies 

E w = Jk 2 + n' k , E a 



1 - JV($, $* ;T,p) — JV($, ; T, /x) 



exp^^/T)-!' 



P + O fc + 2pfi fc: 



where the primes denote derivatives with respect to p = (a 2 + ir 2 )/2 of O = fl + ca. The fermion distribution functions 
AT($,$*;T,/i) and jV($,$*;T,/i), 



M$ $*T = 1 + 2$* exp[/3(£ g - M )] + $ exp[2/3(g g - M )] 

1 ' ' ' W l + 3$exp[2^(^-/i)]+3$*exp[/3(E g -/i)]+exp[3y9(£; ? -/x)] : 
7V($,$*;T, M ) = Ar($*,$;T,-/x), 

are modified because of the coupling to the gluon field. The quark energy is given by 



E q = ^k 2 + 2g 2 p. 



(9) 
(10) 

(11) 



The function q\ is defined as q\{n, q, s) = (2n+l—2s)\q\B and has the same structure as in mean-field approximation, 
see Ref. [5]. 

If one replaces the energies of particles by their tree-level approximation and integrate the flow equation with respect 
to the scale k, then, after integration by parts, the effective thermodynamic potential f2 reduces to the one-loop result 
for the PQM model in the presence of non-zero magnetic fie ld. This also proves validity of the flow equation |8j). 



The minimum of the thermodynamic potential is de- 
termined by the stationarity condition 



da 



da 



(12) 



We solve the flow equation numerically with the 
initial cutoff A = 1.2 GeV (sec additional details in 
Ref. [TTj). The initial conditions for the flow are fixed to 
reproduce the in-vacuum properties (T = p = 0, eB = 
0) : the physical pion mass = 138 MeV, the pion 
decay constant / ff = 93 MeV, the sigma mass m a = 600 
MeV, and the constituent quark mass m q = 300 MeV at 
the scale k — ¥ 0. We treat the symmetry breaking term, 
c = m 2 /^ as an external field. We also neglect the flow of 
the Yukawa coupling <?, which is not expected to be sig- 
nificant for the present studies (see e.g. Refs. [18l H9] V 

The thermodynamic potential (131 does not con- 
tain contributions of thermal modes with momenta 

larger than the cutoff A. To obtain the correct high- The Polyakov loop variables, $ and $*, are then deter- 



temperature behavior of thermodynamic functions we 
need to supplement the FRG potential with the contri- 
bution of the high-momentum states. For this, we follow 
the procedure described in Ref. [T7]: at high k > A the 
meson contribution to the flow in the equation is disre- 
garded and only the flow of massless quarks interacting 
with the Polyakov loop is considered. 

By solving Eq. ([8]), one obtains the thermody- 
namic potential for the quark and meson subsystems, 
Qk^o{$, $*; T, p), as a function of the Polyakov loop 
variables <£> and $*. The full thermodynamic po- 
tential f2($,$*;T, p) in the PQM model, including 
quark, meson and gluon degrees of freedom, is ob- 
tained by adding the effective gluon potential W($,$*) 

to fik_> (*, 

n($,$*;T, M )=fi fc ^ ($,$*;T,M)+W($,$*). (13) 
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mined by the stationarity conditions: 

^O($,<i>*;7» = 0, (14) 
A-n(<S>,<S>*;T,n) = 0. (15) 

B. The mean-field approximation 

To illustrate the importance of meson fluctuations on 
the thermodynamics of the PQM model we compare the 



FRG results with those obtained in the mean-field ap- 
proximation. As we earlier alluded to the one-loop ther- 
modynamic potential, the mean-field approximation can 
be also directly obtained from the flow equation Q by 

(i) neglecting fluctuations of the meson fields and replac- 
ing them by their classical expectation values; 

(ii) integrating the flow equation with respect to the scale 

(iii) performing integration by parts; 

(iv) and, finally, by minimizing the thermodynamical po- 
tential with respect to the mean fields. 



In this way, we will obtain the thermodynamical potential for the PQM model in the mean-field approximation, 
which coincides with the direct mean- field treatment (see Ref. [3]) for the Lagrangian |l]): 



Qmf = **) + U((a),(n) = 0) + Sl q9 ((<r),*, $*). 
Here, the contribution of quarks with the dynamical mass m q — g(a) is given by 



(16) 



-r£ E E 

f=1.2 g=_i I n=0 



\Q ff \B f d Pz f N C E, 



2n 



2tt 



T 



^ + In gM ( (a) , 4>, <T ; T, M ) + In g^ ((<r),$,$*;T,fi) \ , 



where 



(17) 



«?W ((cr), <&,$*;!» = l + 3$exp[-(^-M)/T]+3$*exp[-2(£; 9 - M )/r]+exp[-3(£; g -M)/T], (18) 
9 (-)«a),$,$*;T,/z) = ^+)((a),$*,$;T,- M ) (19) 



and i?q = \/p 2 + q±(n, Qff> s ) + m o is the quark quasi-particle energy (see the definition of q\ in the previous section). 



The first term in Eq. ([17]) is a divergent vacuum fluctuation contribution, which has to be properly regularized. 
Following Rcfs. [11] and [3j, we employ dimensional regularization to obtain: 



2tt 2 



/=1>2 



(20) 



where M is the renormalization scale, 



i 2 q /(2\QffB\) and £(t, x) is the Hurwitz zeta function. 



In Eq. (20) the prime denotes the derivative with re- 
spect to the first argument ('(t,x) — 8Q(t,x)/dt. The 



divergent contribution were subtracted from Eq. ( 20 1 to 
reproduce B — > result for the vacuum contribution [11] 



n vac = 



N r N 



f 4, ( m q\ 



(21) 



Indeed, using the asymptotic expansion of C'( — ^-i x f) m 
the limit B — > or equivalently Xf — > oo (see Ref. [20] ) 



T 2 
X f 



X f 



hixf 



4 



1 

12 



Oix- 1 



f i 
(22) 

it is easy to prove that liniB^o — 

In case of eB — 0, the relevance of the vacuum con- 
tribution for the thermodynamics of chiral models was 
demonstrated and studied in detail in Refs. [TT] and [3T] . 



As usual, the equations of motion for the mean fields 
are obtained by requiring that the thermodynamic po- 
tential is stationary with respect to changes of (a), $ 
and $*: 



MF 



MF 



MF 



8(a) 



9$ 



<9$* 



(23) 



The derivative of the thermodynamic potential with 
respect to the chiral order parameter a also involves 
8Q'[—l,Xf)/8xt which is given by 



dC(-l,x f ) _ 1 
dx f 2 



In 



r(x) 



(24) 



The model parameters are fixed to reproduce the same 
vacuum physics as in the FRG calculation. 
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FIG. 2: The derivative of the order parameter with respect to T as a function of the temperature for the different values of 
the magnetic field eB measured in units of m\. The left (right) panel shows the FRG (mean-field) results. 



III. PHASE DIAGRAM IN THE T — eB PLANE 

In this section, we explore the properties of the chiral 
crossover transition at finite temperature T and magnetic 
field eB in the functional renormalization group (FRG) 
approach and in the mean-field approximation for the 
PQM model. The comparison of the results in two ap- 
proaches allows to pin down the role of meson fluctua- 
tions. As our calculations show meson fluctuations lead 
to even steeper rise of the transition temperature with in- 
creasing field and, thus, are unable to describe the qual- 
itative discrepancy between recent LQCD results [9] and 
the mean-field model predictions, which have been qual- 
itatively confirmed in early LQCD calculations [5] . 

In the FRG approach, the thermodynamic poten- 
tial ( 13 ) at finite temperature and chemical potential is 
obtained by solving the flow equation ^ numerically, us- 
ing the Taylor expansion method (for details see Ref. [17] 
and references therein) . This approach has been success- 
ful applied to the thermodynamics at finite density and 
temperature [TTI [21] . 



The structure of the chiral phase diagram can be re- 
vealed by investigating the chiral order parameter a . In 
Fig. 1, the dependence of the order parameter on the 
temperature normalized by T pc |32j for different values 
of the magnetic field eB is shown in the FRG approach 
and the mean-field approximation. Our results are in 
agreement with those of Refs. [TTJ [25] at zero magnetic 
field and lead to the same conclusion: meson fluctuations 
result in a strong smearing of the temperature depen- 
dence of the order parameter, decreasing the strength 
of the transition. By considering the derivative of the 
order parameter with respect to the temperature, see 
Fig. 2., we conclude that the transition in both models 
shifts to higher temperatures with increasing magnetic 
field. Figure 2 also demonstrates that the strength of the 
transition characterized by the peak in \da/dT\ increases 
slightly with magnetic field. This increase is minute for 
the mean-field approximation and more pronounced if 
meson fluctuations are taken into account by the FRG 
approach. 

The same conclusion on the minor modification of the 
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FIG. 3: The mass of the a as a function of the temperature for the different values of the magnetic field eB measured in units 
of m\. The left (right) panel shows the FRG (mean-field) results. 



transition strength follows from Fig. 3, where the mass 
of c-meson as a function of temperature and magnetic 
field is shown. Figure 3 demonstrates, that in the FRG 
approach, the mass of c-meson in the broken phase is 
modified by the magnetic field stronger than the one in 
the mean-field approach. This is expected because the 
charged meson fluctuations included in FRG make the a 
mass more sensitive to the magnetic field. 

At the chiral critical end point, the mass of er-meson 
(order parameter) is vanishing m a — > 0. Consequently, a 
mass decreases along the crossover line towards the crit- 
ical end point in the T — /i^ plane. Contrary to this at 
zero hb and non-zero magnetic field eB, we do not ob- 
serve this trend. In both approaches, the value of a mass 
at the minimum is almost independent of the magnetic 
field. We checked this fact up to very high magnitudes 
of the magnetic field eB = 30m^. Based on this ob- 
servation, possible chiral critical end point in the PQM 
model can be excluded in a very wide range of the mag- 
netic field [33]. After considering the phase diagram in 
the T — eB plane, we will return to this problem again 
addressing it from a different side. 

By locating the maximum of the peak position in 
\d<j / dT\ as a function of T at a given B we compute 
the phase diagram of the chiral crossover transition in 
the T — eB plane. The phase diagram is shown in Fig. 
4. Contrary to the recent lattice findings [5], the slope 
of the transition line, £, is positive. The inclusion of 
meson fluctuation increases the slope calculated at large 
values of the magnetic field from £mf ~ 4.7 • I0~ 3 to 
£frg ~ 7.1 ■ I0~ 3 . Performing the parametrization 



T pc (B) 
T 

1 nr 



= 1 + A 



eB 

pc 



(25) 



of the phase transition line, which was introduced in 
Ref. [8], we obtain «mf = 6 • I0~ 4 in the mean- field ap- 
proximation and aFRG = 2.7- 10~ 3 in the FRG approach. 
In the PQM model, the role of meson fluctuations at a 
finite magnetic field may be understood by the follow- 



ing considerations. At zero magnetic field, the meson 
contribution to the flow at any temperatures (including 
T = 0) reduces the chiral condensate, i.e. works towards 
the chiral restoration. At a finite magnetic field, the 
charged pions acquire an additional effective mass pro- 
portional to eB, which penalizes meson contribution to 
the restoration of the chiral symmetry at non-zero mag- 
netic field. The quark contribution is not that trivial 
for analytic considerations: while the finite temperature 
part acts towards the restoration of the chiral symmetry, 
the vacuum part has an opposite effect. The numerical 
mean-field calculations show that the quark contribution 
increases the transition temperature. 

Another interesting issue that can be studied in LQCD 
calculations and that may shed light to the QCD phase 
diagram in the three dimensional T—eB—fig space is the 
curvature of the crossover transition line in T — /is plane 
k for non-zero magnetic field. The curvature is defined 
by 



pc(-B)Mb) 



r pc (s,/i S = 0) 



= 1 + k(B) 



Mb y 
T ) 



-o 



Mb 
T 



(26) 



In Fig. 5, we show the curvature of the crossover line 
k(B) normalized by its value at zero magnetic field kq — 
k(B = 0) [34J . The curvature in the FRG approach at 
higher magnetic field shows stronger dependence. 

The issue of the chiral critical end point and subse- 
quent chiral first-order phase transition in the T — eB 
plane can be addressed by a completely different ap- 
proach of the QCD inequalities developed in Refs. [2"3T - 
121)] and recently discussed in Ref. [TD]. Background 
electromagnetic field does not break neither positivity 
of the Dirac operator (V) nor its 75-hermiticity [35] 
(75^75 = £> + ). Thus the results of Ref. [TU] can be 
with no modification extended to the case of the non- 
zero magnetic field. The QCD inequalities 23 26 can 
be translated to those for the meson masses m a > m ff , 
where m T is the mass of the lightest pseudoscalar pion 
and m a is the lightest meson mass in the channel a, 
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consult Ref. [TU] for details. This inequality is rigorous 
for non flavor singlet a. In the large N c limit, however, 
owing to 1/N C suppression of the flavor disconnected di- 
agrams, this inequality becomes applicable for flavor sin- 
glet channel too. For us it is essential, that m a > m n at 
the leading order of large N c expansion. From the last 
inequality, it follows that at any finite value of the pion 
mass the correlation length proportional to l/m a is finite, 
i.e. the chiral second-order phase transition is impossible. 
For the sake of argument, we remind the reader that it 
has been rigorously established in high precision LQCD 
calculations of different groups 27-29] that for the phys- 
ical pion mass and at zero magnetic field and T rs 155 
MeV [29j [30] the chiral transition is a smooth crossover. 
In principle, the chiral crossover transition may turn into 
the first-order one at a finite value of the magnetic field. 
This, however, may happen only via second-order criti- 
cal end point. The above mentioned argument disfavours 
such possibility at least in the large N c limit. 

We note, however, that this argument applies only to 
the chiral transition. Existence of other phase transitions 
(e.g. the first-order deconfincment phase transition) at 
high eB cannot be excluded. 



IV. SUMMARY AND CONCLUSIONS 

Results of the recent LQCD calculations on the depen- 
dence of the transition temperature on the magnetic field 
disagree already at a qualitative level with those obtained 
previously, as well as with various low energy effective 
model of QCD. In this article, we addressed this issue 
in the Polyakov loop-extended quark-meson model be- 
yond the mean-field approximation. We showed that the 
inclusion of meson fluctuations, which presumable were 
suppressed in the early LQCD fluctuations and were ne- 
glected in the mean-field models, is unable to resolve the 
above mentioned disagreement. 

We calculated the phase diagram of the Polyakov loop- 





FIG. 5: The curvature of the transition line in T-/j_b plane, 
ft, as a function of the magnetic field. 



extended quark-meson model in the mean-field approxi- 
mation and in the functional renormalization group ap- 
proach. Both approaches result in a shift of the tran- 
sition temperatures to higher values then that at zero 
magnetic field. Moreover, the relative increase of the 
transition temperature is larger if meson fluctuations are 
taken into account. 

Although we observed that the transition strength in- 
creases with increasing magnetic field, we see no evidence 
in favour of possible chiral first-order phase transition at 
finite values of the magnetic field eB. 

Based on the large N c non-go theorem of Ref. [TUj , we 
provide another indication against a chiral critical end 
point and a change from the chiral crossover to a first- 
order phase transition in the T — eB plane for non-zero 
pion mass. 

Finally, in the PQM model calculations, it was also 
shown that the magnetic field increases the curvature of 
the transition in the T — \lb plane. 

In this model, we do not take into account a contri- 
bution of the charged vector mesons, e.g. p . As was 
shown in Ref. [31], in a high magnetic field (eB c ss m 2 p ) : 
the p-meson condensate may form and drastically change 
properties of nuclear matter. 



FIG. 4: The phase diagram in the FRG approach and mean- 
field approximation. 
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